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Abstract 

The Lipkin-Meshkov-Glick (LMG) model has a Schwinger boson realization in terms of a two-level boson 
pairing Hamiltonian. Through this realization, it has been shown that the LMG model is a particular case 
of the SU (1, 1) Richardson-Gaudin (RG) integrable models. We exploit the exact solvability of the model to 
study the behavior of the spectral parameters (pairons) that completely determine the wave function in the 
different phases, and across the phase transitions. Based on the relation between the Richardson equations 
and the Lame differential equations we develop a method to obtain numerically the pairons. The dynamics 
of pairons in the ground and excited states pairons provides new insights into the first, second and third 
order phase transitions, as well as into the crossings taking place in the LMG spectrum. 

Keywords: Richardson-Gaudin models, Lipkin-Meshkov-Glick model, Quantum phase transitions, 
Integrable models 



1. Introduction 

The LMG model was introduced in nuclear physics to mimic the behavior of closed shell nuclei [TJ. It 
is a simple model with one quantum degree of freedom. The dimension of the Hamiltonian matrix increase 
linearly with the size of the systems allowing its exact diagonalization for large system sizes. As such, the 
model has been extremely useful to test many-body approximations to nuclear problems (see for example 
[2] and references therein). More recently the model found applications to many other areas of physics 
like quantum spin systems [3], ion traps [4], Bose-Einstein condensates in double wells [5] or in cavities 
[B], and in circuit QED [7j. The model has been also utilized to study quantum phase transitions (QPT) 
[51 |S] and their relations with quantum entanglement properties (TUJ HI] , as well as to explore excited states 
QPT [12] and quantum decoherence [T3]. On a different respect, the LMG model was shown to be exactly 
solvable [H] and quantum integrable [Jjj] as a particular limit of the SU(1, I) Richardson-Gaudin integrable 
models [161 117j . The most important feature of the exact solution is that it provides a unique form for the 
wavefunction of the complete set of eigenstates of the model in terms of a set of pair energies or pairons 
obtained as a solution of the non-linear coupled Richardson equations. The distribution of pair energies in 
the energy space change dramatically close to a critical point. A typical example is the exactly solvable 
Px + ip y fermion superfluid derived from the SU{2) hyperbolic RG model [TH] [TH] [2U] . The model has two 
interesting lines in the phase diagram of density versus coupling constant: a) the Moore-Read line in which 
all pairons collapses to zero energy, and b) the Read-Green line in which all pairons are real and negative. 
While in the first case the existence of a QPT is still debated, in the second case it has been shown that 
Read-Green line corresponds to a third order QPT. 

The LMG can be mapped to a two-level boson systems by means of the Schwinger representation of 
the SU(2) algebra. In this representation the LMG Hamiltonian transforms to a two-level boson pairing 
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model associated with a SU (1,1) ® SU (1,1) algebra [H]. Within these models, the relation between the 
distribution of pair energies and the occurrence of a QPT has been discussed in Ref. |21j in connection with 
the s-d dominance in the Interacting Boson Model of nuclear physics. More recently, a thorough analysis 
of critical points of the two-site Bose-Hubbard model in terms of the roots of the Richardson equations has 
been presented |22j . showing the intimate relation between quantum criticality and the rapid change in the 
behavior of the pairons. In this paper we will continue these studies focusing on the generalized LMG model 
which has a rather rich phase diagrams with lines of first and second order QPT and a triple point with a 
third order QPT. Moreover, the pair energies in a region of the parameter space display a behavior similar 
to the p x + ip y superfluid model, opening the possibility of correlating the physics of both models in the 
critical regions. 

We will start by introducing the LMG model, its Schwinger boson representation leading to a two-level 
boson pairing model, and the mean field phase diagram of the model in section [2j In section [3] we will 
introduce the SU(1, 1) RG models and discuss the limits leading to the LMG model. We will introduce in 
this section a robust numerical method to solve the Richardson equations based on their relation with the 
Lame ordinary differential equation. Section [4] is devoted to the study of the behavior of the ground state 
pairons close to the phase transition, and the region in parameters space which shows a behavior similar to 
the Px + ipy model. Finally, in section[5]we will describe the RG solutions for the excited states. Concluding 
remarks are given in section [6] 



2. The Lipkin-Meshkov-Glick model and its bosonic representation 

The LMG model is based on the SU(2) algebra, whose three elements satisfy the commutation relations 

[S +) S-]=23 S , [S Z ,S±] = ±S±. 

The three elements can be considered as the components of the pseudo-spin operator S. They commute 
with the Casimir operator of SU(2), S 2 = ±(S + S- + SS+) + S 2 Z . In terms of these elements the LMG 
Hamiltonian can be written as 

H L = eS z + ^ (S 2 + + S 2 _ ) + | (S+S- + SS+) , (1) 

Note that H L does not commute with the z component S z for A ^ but it commutes with the total pseudo- 
spin Casimir operator S 2 . Therefore, the Hilbert space of the model can be separated in different sub-spaces 
labeled by the eigenvalues of the total pseudo-spin j(j + 1), with basis Hj = {\jm) : m = —j, —j + l,...,j — 
1, j}. Additionally, the LMG Hamiltonian commutes with the parity operator P = expi%(S z + j), yielding, 
for a given j, two invariant sub- spaces (P = + and P = — ), which are spanned, respectively, by the basis 
Hj+ = {\jm) : m = —j, —j + 2, — j + 4, ...} and %j- = {\jm) : m = — j + 1, — j + 3, — j + 5, ...}. From now 
on and for the sake of simplicity, we will assume integer values for j. The semi- integer case can be worked 
out following the same lines with some slight modifications. For integer j the dimensions of the invariant 
subspaces, 1-Lj+ and are j + 1 and j respectively. 

Having introduced the LMG Hamiltonian in terms of SU(2) operators, a physical realization of the model 
requires a representation of the algebra either in terms of a collection of spins or in terms of a fermionic 
or bosonic system. In its original presentation the SU (2) operators were expressed in terms of a collection 
of 2j fermions distributed in two levels, each having a 2j fold degeneracy. Instead we will make use of 
the Schwinger boson representation of the SU(2) which allows to a simple connection with the bosonic RG 
integrable models. 

The Schwinger representation of the SU(2) algebra in terms of two bosons is 

&tft — a +a n b -h a f f 

S z = = , S + = b J a, S-=a , b, (2) 

with a and b boson operators, satisfying the usual commutation rules [a, <v] = [b, w\ = 1 and [a, b] = [a, b'} = 
0. Inserting the boson mapping ([2| into the Hamiltonian ([!]) the bosonic version of the LMG Hamiltonian 
reads: 
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H L = V^^b + ^a^a + ± (tf 2 a 2 + a^b 2 ) + 7 (bWab) 



(3) 



Using the Schwinger representation, a basis for the Hilbert space with total pseudo-spin j can be written 



in terms of boson creation operators as \jm) — \n a = j — m, n b = j + m), where \n a , n b ) 



(a t )"°(fe t )" b 



|o>, 



\/n a !iij,! 

with |0) the boson vacuum. Note that for a given j the total number of bosons is constant N = n b + n a = 2j. 
Likewise, the positive and negative parity basis in the Schwinger representation are given byin" 

K j+ = {\n a = 2j, nb = 0), \n a = 2j - 2, n b = 2), \n a = 2j - 4, n b = 4), ...} 

Hj- = {\n a = 2j - l,n b = 1), \n a = 2j - 3,n b = 3), \n a = 2j - 5,n b = 5), ...}. 

A detailed analysis of the LMG phase diagram in terms SU{2) coherent states with definite parity has 
been performed in Ref. [5j. The different phases of the model and the order of their transitions were 
identified. Here we repeat that analysis using the Schwinger representation and the boson coherent state 



|0), 



(4) 



where z are c- numbers parametrized as z = pe 10 . The expectation value of the LMG Hamiltonian ([3|) in the 
coherent state \z a z b ) = is 



(z b z a \H L \z a z b ) 



7' 
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^ + (Acos(2(^-^))+ 7 KP6- 



(5) 



The constraint n a +n b = 2j implies that coherent states parameters should fulfilled p 2 + p 2 — 2j. Enforcing 
this relation, the energy surface is 



E[p b ,9] = (z b z a \H L \z a z b ) = — 



2je 



2J-1 



A 



2j-l+jBt 



jB 
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with A 



7x+7i,-2(2j-l) 



re-scaled parameters defined in [5] 



and Be = + 7^ + (7^ — 7l/ ) cos 26*^ , where we have used 9 = 6 a — 9 b , and 



(7*, 7y) = — — -(7 + A,7-A). 



(6) 



the 
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In the thermodynamic limit, j — > 00, the energy per particle, £ [p b ,0] = E[p b ,9]/(2j), simplifies to 



2£\p b ,( 



1=2 + 5, 



Bo 



(8) 



where terms of order OJl/j) have been neglected. The phase diagram of the LMG model is obtained by 
minimizing the energy (8 1 with respect to the variables 9 6 (— 7r,7r] and |j € [0, 1] for different values of the 
model parameters, Jx-Jy The different phases, separated by dashed lines in FigJT] are described in Table I 
where we have classified the values of the parameters 9 and p b characterizing the coherent state Q at the 
absolute minimum. Additionally, Table I shows the energy per particle and the expectation values of the 
operators S 2 and S 2 that play the role of order parameters. The critical line j x = j y < — 1 is a special case 
because the relative phase 9 is completely undetermined, i.e. the minimum in energy is independent of 9. 

We are now ready to establish the phase diagram of the LMG model in the Schwinger boson representa- 
tion. In complete accord with Ref.[8 , three phases are identified, which are distinguished by the occupation 
of boson b and the relative phase (9) of the coherent state parameters z a and z b . Alternatively we can char- 
acterize the three phases by the order parameters ^|f^ = (p 2 /j)(p 2 /j) cos 2 9 and < jf-) = (p 2 /j)(p 2 /j) sin 2 9, 
where we have neglected terms of order 0(1/ j). Phase A (j x > —1 and 7 y > —1) has zero occupation of the 
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Table 1: Phases and their order parameters in the LMG model 



boson b, (f^J = ^ ZaZb \ n ^\ ZaZ ^i — o. Therefore the two order parameters are also zero, (S x ) = (S%) = 0. 

\ •* J min 

In phase B, (7,, > 7^ and ~f x < —1), the coherent state mixes a and b and the order parameter (S x )/(j 2 ) is 
finite. Finally, phase C with (7^ < j x and 7^ < — 1) is the mirror of phase B corresponding to an exchange 
between x and y. Upon inspection of the order parameters, we can immediately recognize that the tran- 
sitions between phases (A-B) and (A-C) are continuous in the order parameters, defining a second order 
phase transition. The transition between B and C is discontinuous in the order parameters characterizing 
a first order phase transition. These facts were confirmed in Ref. [H] by analyzing the energy derivatives. 
At the triple point, 7^ = 7j, = —1, both order parameters converge to avoiding the discontinuity of the 
first order critical line. As shown in [8] this critical point represents a third order phase transition when it 
is traversed in the direction indicated by the arrow in Figjl] 



3. SU(1, 1) Richardson-Gaudin and LMG models 

While the exact solution of the LMG was derived in Ref. [T3] using an algebraic approach based on the 
Bethe ansatz, the connection between the LMG Hamiltonian and the SU(1, 1) RG models was established 
later [15j . In order to make the present work self-contained and to fix the notation, we will here derive the 
exact solution of the LMG Hamiltonian from the more general SU (1,1) RG models following a different 
path. 

The non-compact SU(1, 1) algebra, defined in terms of ladder (K+ , K_ ) and weight (K z ) operators 
resembles that of the SU(2) group, differing in a sign in the commutation relations 

[K^K^] = ±K ± , [K+,K-] = -2K Z . 

Let us now consider N c different copies of the SU(1, 1) algebra, and construct N c linear and quadratic 
hermitian combinations of the three elements of the algebra 



i? 4 = ^-2 5 ]T 



X(t i: tj) 



(K+KJ+KTK+) - Z(t h tj)KtK] 



(9) 



where i, j label each of the N c copies and g is an arbitrary parameter. The structure of the operators Ri 
is such that they commute with the total K z operator (K z = J^i-K-i)- ^ has be shown [T5] that the N c 
operators commute among themselves ([Ri, Rj] = 0), defining and integrable model, if the functions X(ti, tj) 
and Z(ti,tj) are anti-symmetric functions of an arbitrary set of parameters U 



V /(1 + ,<)(1 + ^) 



1 + std 



X(U,t j ) = l , Z(t l ,t J ) = - (10) 



"31- f-t- ' t -t 



The parameter s can take three different values s = 0, 1, —1, defining the rational, the trigonometric, and 
the hyperbolic families of SU(1, 1) RG models respectively. 
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Figure 1: Phase diagram of the LMG model and Richardson-Gaudin areas of integrability in the 7^-7y parameter space. The 
light quadrants correspond to hyperbolic model (s = —1) with light gray for < t < 1 and white for t > 1. Dark gray 
quadrants correspond to the trigonometric model. Upper quadrants correspond to positive g while lower quadrants to negative 
g. The dashed lines separate the three different phases (A,B and C) of the LMG model discussed in the text. The triple point 
(— 1, —1) in the intersection of the lines is a third order transitions when it is traversed in the direction indicated by the arrow. 
The horizontal line -y y = and the vertical one -y x = correspond to the rational version of the model (s = 0). 



The LMG model is obtained in the limit of two SU(1, 1) copies, which we will label as a and 6. In the 
pair boson representation of the SU(1, 1) algebra, the elements of the two copies are 



K+ = ±<Ja1 K~,= \aa, K' a = \foa+^\ 



Kt = \b^ K-,= lbb, K^ 1 -(blb+ 1 -y (11) 

The irreducible representations (irreps) of the non-compact SU(1, 1) algebra are dimensionally infinite, but 
they possess a minimum weight state defined by K~\MW) = 0. For the previous bosonic representation, 
these states are given by \vi = 0) = |0)j and \v.i = 1) = |l)i, where |0)i is the vacuum of bosons i = a,b. The 
parameters i/i are the so-called seniorities of each of the SU(1, 1) copies. The seniority quantum number, 
Ui, counts the number of unpaired bosons i and can take only two values or 1. If = the number of 
bosons i (n^) is even, and odd if v, L = 1. 

Inserting the boson pair representation of the two copies ( 11 ) in the integrals of motion ^ we construct 



the two integrals of motion of the LMG model. We can verify the the sum of both integrals gives the 
conserved quantity K z . Taking the difference between both we obtain: 

R b - R a = Kl - Kl - 2gX ba [K+K~ + K+R-] + AgZ ba K§K z a 

= \ (tfb at a) - (b^a 2 + a^b 2 ) + gZ ba fob + foa + , 

with Z ba = Z(t b ,t a ) and X ba = X(t b ,t a ). Comparing with the LMG Hamiltonian in the Schwinger repre- 
sentation ([3]), one finds the following relation between the LMG model and the integrals of motion of the 
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transition 


line 


RG parameters 


first-order 


1x 7y •> 7a: ^ 1 


s = -l i = l (2j-l)# = -l 


second-order 


lx = -1, 7y > -1 


« = +1 (2j - 1)5 = +7 
s = -l £>1 (2j-l) 9 = -| 


second-order 


lx > -1, ly = -1 


s = +1 (2j - l).g = -t 
s = -l t<l (2j-l)g = -t 



Table 2: Transition lines in phase diagram -y x , jy and their translation to the RG parameters. 



SU(1, 1) RG models 

7 A T 

H L = e(R b - R a ) - with gX ba = --, and gZ ba = -. (12) 

Without any loss of generality, we choose the parameters entering in and i?5 Q as = — i a = i, with 

t > 0. The functions X(, a and Z;, a reduce to X ba — >ft an< ^ = ■ Then, the relation between the 

LMG Hamiltonian parameters (A, 7, e) and those of the SU(1, 1) RG models are 

A 1 + st 2 7 1 - st 2 



e " 9 2t ' e 5 " 2t 



Or in terms of the 7^ and 7 a parameters ^ we have 



(7x,7v)= 2l 7 ^(7 + A,7-A) = (2j-l). 9 (^-^,^ . (13) 

The relation between the j x and j y parameters and those of the RG model classify the quadrants of the 
phase diagram of Figjljin terms of the hyperbolic (s = —1) and trigonometric (s = 1) RG models. The first 
(s = — 1, g > 0) and third (s = — 1, g < 0) quadrants correspond to the hyperbolic RG model, whereas the 
second (s = 1, g > 0) and fourth (s = 1, g < 0) are associated with the trigonometric model. These regions 
are indicated in Figure [lj by dark gray zones for the trigonometric model and light zones for the hyperbolic 
one. The rational RG model is limited to the j x = and 7 y = lines. 

The phase transition lines, discussed in the section [2] and shown in Fig. [l]by dashed lines, are translated 
to the RG parameters in table [2] 

The LMG model has symmetries that relates the spectrum of systems in two different points in the 
parameter space. The first of these symmetries is a point reflection trough the origin, and relates systems 
obtained from a simple change of sign in the parameters (7^,7^) — > ilxily) = ( — 7x> — ly)- This change 
in sign is equivalent to a global sign change in the Hamiltonian, implying that the spectrum of a system 
is minus the spectrum of the transformed system. In terms of the RG parameter, this transformation 
corresponds to g —¥ g' — —g. A second symmetry of the LMG model is a reflection across the line "f y = j x . 
Two mirror points of the phase diagram symmetrically located around this line have exactly the same 
energy spectrum, as a result of the invariance of the SU(2) algebra under the canonical transformation 
S+ —} —iS+, S- —t iS-, S z — > S z , corresponding to b — > ib, a — > a in the Schwinger boson realization [23 . 
In terms of the RG parameters, this symmetry implies that systems with parameters related by 

(<7,t)^ (</,*') = (-«<?,!/*), (14) 

have the same spectrum. 

The 7 y = j x line (s = — 1 and t = 1 in terms of RG parameters), located within the hyperbolic regions of 
the phase diagram, has the peculiarity that the RG solutions display a singular behavior as it will be shown 
in section |4| However, the eigenstates along this line can be easily obtained in closed form by resorting, for 
instance, to the LMG Hamiltonian in terms of pseudospin operators. The condition j x — j y implies A = 



G 



and results in the following LMG Hamiltonian 

Hl = £^S Z + - 2 J C _ - (S+S- + SS+^j =e(S z + g (S 2 - 5 7 2 )) , (15) 
which commutes with both S 2 and S z , and has eigenstates \jm) and eigenvalues 

E m = e (m + 7 * (j(j + 1) - m 2 )\ = e(m + g (j(j + 1) - to 2 )) . (16) 

The conservation of the S z operator along this critical line allows the existence of real crossings between states 
of the same parity. These crossings take place when E m = E m ±2- For instance, the P+ ground state energy 
for g = 0, E m= ^j, crosses the first P+ excited state energy, i? m= _ J+2 , when -f x — j y = — (2j — l)/(2j — 2), 
or in terms of RG model parameter when s = —1, t = 1 and g = —l/(2j — 2). 

Having related the two-level SU(1, 1) RG models with the LMG model, we can now explore the exact 
RG solutions in each subspacc defined by the number of boson pairs M and seniorities v a , v b- For integer j 
the number of Schwinger bosons, N = n a + rib = 2j, is even, the seniorities are equal, v a = Vb = v = or 1, 
and the number of boson pairs is M = j — v. The seniority sectors v = and v — 1 correspond to the two 
invariant sub-spaces P = + and P = — respectively 

The eigenvalues r^, of the Ri integrals of motion are [15] 

(M 
1 + 2 Z (U, e a ) + 2^2 djZtjk, tj 



where di — (1/2) (^ + |), Z is the function defined in (10 1, and e a are the so-called spectral parameters or 
pairons. Each particular eigenstate is completely defined by a particular M pairon solution of the coupled 
set of Richardson equations 

M 

1 + 2 d t Z{e a , U) + 2^2 Z(e a ,e p ) = 0. 

For the particular case of the LMG model, with two SU (1,1) copies, tb — —t a — t, and integer j 
(i> a = Vb = v = 0, 1), the Richardson equations reduce to 

1 + se 2 + 9 \ + 2~) \e a + t + e a - t) + 2 ' 9 ^ e a - ep °' 



The eigenvalues of the LMG Hamiltonian are given by 

E L = <r b -r a )-l=ge^-^^ + l) + 2 g eL+l)tj: 1 -^§. (18) 



4 3 2t V 2 J ^ t 2 

The unnormalized eigenvectors common to the two integrals of motion (Ri) and, consequently, to the LMG 
Hamiltonian are 



a=l 



The Richardson equations can be interpreted as an electrostatic problem in two-dimensions [2H US] . In 
order to make explicitly this connection we rewrite the Richardson equations as 

Qc Qd , v + 1 ( 1 , 1 \ , ^ 1 



P c e a -P D 2 \e a +t e a -tj e a - ep 



E — - = 0. ( 2 °) 
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with the effective charges Qc and Qd 



(1 1 2j ~ 1 

Qc — 



4g^-s 4 



Qd = --I-**-!, (21) 



located at position Pc = —l/y/^s and Pd = l/y/—s respectively. The pairons have a positive unit charge 



and they are located at positions e a in the complex plane. Eq.(20) describes the electrostatic interaction 



of a set of M pairons with positive unite charge in a two dimensional space. The first two terms in (20 1 



describe the electrostatic interaction of the pairons with the two charges Qc and Qd- For the trigonometric 
case (s — 1) the effective charges are complex Qc = — 2 ~' 4 ~ 1 — J^, Qd = Qci an d they are located in 
Pc = i, Pd = —i- Whereas in the hyperbolic case (s = —1) both, the charges and their positions, are real, 
Qc = — 4- 1 + and Qd = — 2j ' 4 ~ 1 — j^, located in Pc = — 1 and Pd = 1 respectively. The third term 



1/4- i 

in (20) represents the interaction of the pairons with two charges at positions =ft. Finally, the fourth 



term corresponds to the mutual repulsion between pairons. Each independent solution of the Richardson 
equations determines the equilibrium position of the pairons in the complex plane. The electrostatic mapping 
will be useful to interpret the pairons distribution in each of the quantum phases of the LMG model and 
the structural changes that take place close to the quantum phase transitions. 

3.1. The Richardson solution as the roots of a generalized Heine- Stieltjes polynomial 

The standard way to solve the Richardson equations is to start from the weak coupling limit where the 



solution is known [see Eq. ( 33 ) below] . The coupling strength g is increased gradually, using the previous 
solution as an initial guess to solve the equations for the updated g by means of a standard Newton-Raphson 
method. A recursive use of this strategy allows to reach the solution for an arbitrary value of g provided 
one is able to develop a method to treat the numerical instabilities appearing when two or more pairons 
converge at the same point (at the position of the U parameters in this case) generating singularities in the 
equations [26 . Recently, two related methods for solving the Richardson equations have been presented 
( [271 12"5"] and [29 ). Both methods exploit the relation between the Richardson equations and the Lame's 
Ordinary Differential Equation which has a generalized form of the Heine-Stieltjes polynomials as a solution. 
The roots of these polynomials are precisely the spectral parameters or pairons. These methods have the 
advantage of being numerically more stable for system of moderate sizes. For larger systems they have 
instabilities due to the large precision needed to calculate the roots of a polynomial of high degree. Here we 
follow the method of references [37J HH] , which is more adequate for systems with a small number of levels 
|30j . as it is the case of the LMG model. 



We begin with the Richardson equations in the form (201, which can be written as 



E— ^ = -E^^> ( 22 ) 

with ( Pk ,Vk) = (Qc, -l/V 1 ^), (Q D , ((2v + l)/4, -t), and (O + l)/4, t) for k = 1, 2, 3, 4 respec- 

tively. 

Let us now define the polynomial P(x) = Yi^=i( x ~ e a) which can be expanded in powers of x as 

POr) = J> fc x fc , (23) 

k 

whose roots are the set {e Q } for a particular solution of the Richardson equations. This polynomial is 
a generalized Heine-Stieltjes polynomial that satisfies the following Lame's ordinary differential equation 



( Appendix A ) 



A(x)P"(x) + B(x)P'(x) - V(x)P(x) = 0, (24) 
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where the functions A, B and V are polynomials defined as 



A(x) = U(o 



Ik) 



B{x)=A{x)Y, 



%Pk 



k=l 



T]k 



with 



A(x) 



P{x) 



and V(x) 



E 



i l^k 



(25) 



(26) 



The polynomials A(x) and B{x) of degree 4 and 3 respectively, depend only on the parameters of the LMG 



Hamiltonian (13). V(x), the so called Van Vleck's polynomial, is at most of third order and depends on the 



values A(77i), which in turn depend on the set of pairons e a , 



V(x)=J2 b > xt - 



(27) 



For a general problem [27] , one can insert the polynomials V(x) ( 27 ) and P(x) ( 23 ) in the ordinary differential 
equation ( 24 ) and, by equating to zero the coefficients at each order in x, one obtains two systems of equations 



for the coefficients bi and . The first set of equations is linear allowing a solution in which coefficients bi are 
expressed in terms of the a^, leaving a second set of non-linear equations for the cii coefficients. Finally, the ai 
coefficients determine the polynomial ( |23[ ) whose roots are the pairons e Q of the Richardson equations. For 
the particular case of the LMG model the first system of linear equations allows to determine the coefficients 
of the Van Vleck polynomial bi (except b ) directly in terms of the parameters of the problem (i.e. they are 
independent of the ai parameters). As a consequence, the second set of equations is linear in the coefficients 
aj. From Eq. (25 I we obtain: 



A{x) 
B(x) 



-sr + (s - t 2 )x 2 + x 



-st 



+ \t 2 (2j - 1) + s(2v + 1)1 x + -x 2 - 2(M - l)x 3 
9 9 



(28) 



After substitution of these polynomials into the differential equation ( 24 ), and from the terms of order M+3, 
M + 2 and M + 1, we obtain 6 3 = 0, 62 = —M(M — 1) and 6 X = sM / g, i.e the Van Vleck polynomial is 
completely determined except for the order zero coefficient: 

V(x) = 6 + —x ~ M(M - l)x 2 . 
9 

We can derive the b Q coefficient and the parameters a, from the orders to M of the differential equation 



(24 1. The result is an eigenvalue equation: 



a 1 



^ D kk >a k , = b a k with k = 0, 1, M, 



k'=0 

where the matrix D kk / is completely determined by the parameters of the model (t,g, M, v, s). Its non-zero 
matrix elements are given by 

D kk _ 2 = (fc-2)(fc-l-2Af)+M(M-l) 

D kk - t = s(k-M-l)/g 

Dkk = H(2j - k)t 2 + s{2v + k)) 

Dkk+i = s(k+l)t 2 /g 

D kk+1 = -s(k + 2)(k + l)t 2 . 



(29) 
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Figure 2: Ground state pairons as a function of g(2j — 1) for the trigonometric case (s = 1) with t = 1 and j = 15, corresponding 
to ■jy = —■yx as indicated with the arrow line in the inset. Here the arrow corresponds to increasing values of g. The vertical 
dashed line indicates the critical value of g = g cr = —l/(2j — 1). Pairons for positive and negative values of g are related by 
e a — > l/e a . 



Therefore, the coefficients at of the polynomial P{x) are the elements of each eigenvector of the matrix D. 
Once the coefficients are known, the e a roots of the Richardson equations are obtained by finding the 
roots of the polynomial P(x). Each of the M + l eigenvectors of matrix D, defines a polynomial whose roots 
e Q correspond to a particular eigenstate of the LMG Hamiltonian. 

It is important to remark here that the drawback or bottleneck of the method resides in the last step. As 
it is well known, finding the roots of high degree polynomials requires a high precision in the determination 
of the coefficients. Therefore, the number of pairs M is limited to = 10 2 — 10 3 . Conversely, the method 
allows to find directly the pairon roots without resorting to the iterative method of increasing gradually the 
coupling strength with the burden of having to deal with the singularities of the Richardson equations. 

4. Numerical Results for the ground state 

In this section we will present and discuss the numerical solution of the Richardson equations for the 
ground state in the different phases using the trigonometric and the hyperbolic model as required for each 
particular phase. 

Trigonometric quadrants 

It order to illustrate typical results for the trigonometric regions (s = 1), we will consider the line 
with t = 1 as a function of g. This line corresponds to 7 = 0, which cancels the third term in the LMG 
Hamiltonian ([I]). The resulting Hamiltonian is the most frequently used in the literature, also known as 
the Lipkin Hamiltonian. In terms of the scaled parameters this line corresponds to 7^ = —~/ x = (2j — l)g. 
For increasing values of g, we move along the line from the fourth to the second quadrant in the phase 
diagram. According to this diagram, a second order phase transition takes place in the thermodynamic 
limit for j y = ±1, or equivalently for g = g cr = ±l/(2j — 1) (see table [2]). The e a pairons for the ground 
state as a function of the ratio g/g C r = ff(2j — 1) are shown in figure 2, for a system with j = 15. 

As it can be seen in the figure, for g ~ all the pairons are located close to t = — 1. As the strength of g 
is increased they expand in the real axis. For negative g, the pairons are constrained to the interval [— t, t], 
and they behave in a very similar way to the rational case already discussed in the context of the IBM-model 
|21) . The second order phase transition can be interpreted as a localization-delocalization transition. The 
pairons initially localized close to t = — 1, expand to the entire interval [— t,t] in the transition point. For 
positive g the solutions are, except for a sign in the wave function, entirely equivalent to the negative g 
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Figure 3: Real part of the ground state pairons as a function of g(2j — 1) for the hyperbolic case (s = —1) with 4 = 1/2 and 
j = 15, corresponding to -y y = A^y x as indicated with the arrow line in the inset. The arrow corresponds to increasing values of 
g. The vertical dashed line indicates the critical value of g = g cr = — l/(4j — 2) corresponding to the line ■jy = — 1. For positive 
g, successive collapses of pairons in the position Pq = — 1 (horizontal dashed line) of the effective charge Q c can be seen. 



case. As it was discussed in section [3j the Richardson solutions for two mirror points symmetrically located 
around the "f x = j y lines [ (g, t) and (—sg, 1/i) in terms of RG parameters] have the same spectrum and the 
wave functions are related by a canonical transformation b —¥ ib. This symmetry is reflected by a simple 
relation between the pairons with negative and positive g given by e a —> l/e a . It is straightforward to 



show that the energy in Eq.(18l is invariant under the transformation (g, t, e a ) — > (—sg, 1/t, l/e a ), and that 



this transformation produces a change in the relative sign of the two terms appearing in the product wave 



function (19), in agreement with the canonical transformation b — > ib. 

In the trigonometric quadrants, the dynamics of the pairons as a function of the control parameters 
[lx, ly\ take place entirely in the real axis and it is very much like the already known dynamics of pairons 
in the rational boson pairing models [21] . 

4-2. Hyperbolic quadrants 

The hyperbolic regions (s = —1) of the phase diagram offer much richer structures than the trigonometric 
regions. In order to illustrate this issue, we will study a system with j = 15 and t = 1/2, which corresponds 
to the line j y = Aj x in phase diagram of FigjTJ The line traverses the third and first quadrants from below for 
increasing values of g, and has a critical point of a second order phase transition for "f y = — 1, corresponding 
to g(2j -l) = -t= -1/2 (see table [2]). 

Fig. [3] shows the ground state pairon roots as function of g(2j — 1). Similarly to the trigonometric 
results, the pairons converge to —t for g — > 0. For negative values of g the pairons are constrained to the 
interval [t, —t] with the phase transitions (vertical dahsed line) signaled by the derealization of the pairons 
in this interval. For the g-positive case (where no phase transtition is expected) an interesting behavior 
of the pairons takes place. As the coupling g is increased the pairons collapse successively to the position 



(Pc = —1) of the effective charge Qc- In Appendix B it is shown that a necessary condition to have Nq 
pairons collapsing to Pc = — 1 is 

ff = ^2j + i-2J\r c ' (30) 

with N c = 1, M and < g c Nc=1 < g c Nc=2 < ••■ < 9n c =m- 

According to this expression, the first collapse occurs for Nc = 1 (one collapsing pairon) at g(2j — 1) = 
g^(2j— 1) = 1, then Nc — 2 pairons converge to P c = — 1 at g(2j— 1) = g^ftj — l) = (2j — l)/(2j — 3), and so 
on. After the collapse of an even number of pairons a new complex conjugated pair of pairons is created. In 
figure [4j the real and imaginary parts of the pairons are shown for positive g values. The successive collapses 
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Figure 4: The collapses of the ground state pairons as a function of g(2j — 1) for the hyperbolic case of figure [3](i = 1/2 and 
j = 15). Left and right panel show the real and imaginary part of the pairons respectively. The collapses take place in the 
position Pq = — 1 (horizontal dashed line in left panel) of the effective charge Qq at values of g given in | |30| . The dashed 
vertical lines indicate the number of pairons Nc = 1, ...,8 involved in the collapse. After the collapse of an even number of 
pairons a complex conjugated pair is created as can be seen in the right panel of the figure. 



and creation of complex conjugated pairs can be clearly seen. This behavior was completely unexpected in 
bosonic RG models, where pairons have always been constrained to the real axis. 

A particular situation occurs when all pairons collapse to Pc = — 1 for g = g € N = 1. At this particular 



point the exact ground state eigenstate ( 19 1 takes the simple form 



l^=lffi-tTi/^ 

which would be the boson version of the Moore-Read state found for the p x + ip y model jTSJ HH1 HO] > derived 
from the fermionic hyperbolic RG model. Likewise, as in the fermionic model, the energy of this state is 
E = 0. However, within the LMG model it can be shown by exact diagonalization of very large systems 
that the collapse of all pairons to Pc = — 1 is not associated with a ground state phase transition. A subject 
still under debate for the p x + ip y fermionic model jTHl HH1 ED] . 



It is worth mentioning here that the condition (301 of pairons converging to the value Pc = — 1 applies 



equally to the excited states, and that similar collapses in the position (Po = +1) of the effective Qd charge 
occurs for excited states in the g < interval for values given by 

9 = 9 c Nn =~ 2j + 1 1 _ 2ND , (31) 

with Njj = 1, ...,M and g% D = M < 9n d =(m-x) < — < 9n d =i < 0- These issues will be discussed in section 
[5] where it will be shown that, even if the collapses are not associated to a ground state phase transition, 
they are related to the crossings of excited states of different parities. 

4-3. The triple point = j x = — 1 

The triple point in the phase diagram of the Lipkin model, (7^, 7j,) = (—1,-1), constitutes one of the rare 
example of a third order phase transition in quantum many-body systems. As such, it deserves a thorough 
study because it could shed light into other third order QPT like the one taking place in the p x + ip y model 
|20j . The third order character of this phase transition reported in [8] is observed when the critical point 
is crossed, for instance, along the line 7 y ~ — j x — 2. In figure [5ja the behavior of the ground state pairons 
close to the triple point is examined for a system of size j — 10, moving in the phase diagram along the 

lines 7 a = — j x + b, for three values of b (b = —2.0 grey solid line, b cr = —2. 1 dashed line, and b = —2.32 
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Figure 5: Pairons close to the triple point (7i,7y) = (—1,-1) for a j = 10 system as a function of t. t = 1 corresponds 
to the line 7^ = ■yy in the phase diagram. The different lines are associated to three values of values of b along the line 
■y y = —■y x + b. The values b = —2.0, b cr fa 2.11 and 2.32 are represented by gray, dashed, and black lines respectively, with 
b C r = — 2(2j — l)/(2j — 2). Panel (a) describes the behavior of the ten pairons in the interval 0.5 < t < 1.0. Panel (b) is a 
close up of the tenth pairons around t = 1.0 showing the discontinuous jump for b = b cr . 



black solid line). We move along these lines using the parameter t of the RG model ( pL3[ ) - The value t = 1 
corresponds to the point 7^=7^ = b/2. The critical value of b is b cr = — 2(2j — 1)/ (2j — 2), which is b cr < —2 
for finite systems and b cr — > —2 in the thermodynamic limit. As already discussed, at the points 7^ = 7^ 
the LMG Hamiltonian is diagonal in the basis \jm = (—j + 2k)) cx (aV) <J ~ fc ) (tftf) k \0) , with eigenvalues 



given by ( 16 1 . The energies Ej m= _j and £ , jm =(-j+2) cross at r y x = j y = — (2j — 1)/ (2j — 2). Therefore, the 
line 7 y = — j x + b cr traverses a point at which the positive parity states \jm = —j) and \j, m = —j + 2) are 
degenerated. 

Fig. [5]a shows the behavior of the 10 pairons in the interval 0.5 < t < 1.0 for the three values of b. In the 
three cases the behavior of the nine lowest pairons is similar, all of them converging to —1 in the limit t — > 1. 
However, the last pairon close to t = 1 distinguishes clearly the three cases studied. For b > b cr the last 
pairon converges to —1 like the other nine pairons, whereas for b < b cr the last pairon converges to +1. The 

critical b — b cr separates both regions. In this case the last pairon converges to e_ = v^2iO=L-^ « 0.865 for 

t = 1. The second panel, Fig. [5|b, shows more clearly the behavior of the last pairon near t = 1 for b = b cr 
and b < b cr . Here, the horizontal scale has been extended to t > 1 using the symmetry transformation 



(14 1 [(g, t, e a ) — > (<?, 1/t, l/e a )]. While for b < b cr the last pairon changes continuously across the value of 
t = 1, for b = b cr a discontinuity in the last pairon at t = 1 occurs due to the crossing of positive parity 
states, jumping from e_ « 0.865 to e + = (l/e_) w 1.156. The degeneracy at the point ^ x — ly — b cr /2 is 
associated to the states \jm = —j) oc (aW)- 7 10) and \jm = (—j + 2)) oc (a^a^)^ _1 ^(&^6^)|0). As it can be 



inferred from the exact wave function) 19 1, the limit t — > 1 with b = b cr produces the right eigenstate 



, , / + twj w( a'a' b^b^ \ . . \jm = — 7) — \jm = (— j + 21) 
|*) = ( a t a t)0-D [—- - + ——J |0) cx |*_> = H — ^, 

which is a linear combination of the two degenerated states. The limit t — > 1 + (with the last pairon converging 
to e + ~ 1.156) produces the other degenerated state, \^+) = (1/V2){\jm = —j) + \jm = (—j + 2))), 
orthogonal to the previous one. In summary, when the system traverses the line 7 y = j x along the line 
j y — —j x + b cr the ground state wave function presents a discontinuity, changing from |^_) to |^+). 

In general, the first order phase transition along the line j y = j x (with 7^ < 1) is due to the crossing of 
the states \ jm) and |j(m + 2)). As it was illustrated above for the particular case of the states \jm — —j) 
and \jm = (— j + 2)), the behavior of the pairons near this line reflects these crossings by a discontinuity 
in their values. As a result, crossing the line j y = j x (with j x < 1) implies a jump from a IvP-m) ground 
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Figure 6: Pairons for the k = 10-th positive parity excited state of a j = 15 system. The perturbative results \33\ are shown 
by dashed lines and the exact ones with solid lines. The total number of pairons is M = j = 15. The k = 10-th state is 
characterized by having 5 pairons close to t a = —t and 10 pairons close to % = +t at weak coupling. Note that the lower 
pairons for negative g are inside the interval [— t,t], and the upper ones outside. The opposite happens for positive g where the 
lower pairons are outside the interval, and the upper pairons are inside. 

state to a l^+m) ground state, where |\t± m ) = ~^{\j m ) ± \j( m + 2))). This first order phase transition 
is signaled by a discontinuous change in the order parameters (S%) and (Si). It can be shown that in the 
thermodynamic limit the order parameters for these states are 



Therefore, a jump from l^-m) to |^ +TO ), for m ^ — j, produces a discontinuity in the order parameters 
characterizing a first order phase transition, in complete accord with the analysis of section 1. For the 
particular case of figure[5](m = —j) both order parameters vanish at the critical point preventing a first order 
phase transition. The critical value for this continuous phase transition is b cr = —2(2j — l)/(2j — 2) — > — 2 
in the thermodynamic limit, corresponding to the triple point j x = 'jy — — 1, in complete agreement with 
the thermodynamic results of reference [S] 

5. Excited states in the hyperbolic LMG model 

We will study in this section the excited states for the hyperbolic regions (s = — 1) of the LMG model in 
terms of the pairon dynamics of the RG model. A similar description for the rational bosonic RG model was 
performed in [21] . Pairons for the rational as well as for trigonometric bosonic RG models are always real. 
The hyperbolic model has the particular feature that the pairons can take complex values. For the sake of 
clarity, let us assume the specific value t — 1/2, with generic results for the cases t < 1. The effective charges 
Qc and Qd of the Richardson equations ( p0| located at positions Pc = — 1 and Pp = 1, are outside of the 
interval [— t, t]. The cases with t > 1, can be inferred from those with t < 1 trough the mirror transformation 



Let us first analyze the limit g — 0. In this limit the LMG Hamiltonian reduces to an one-body 
Hamiltonian H — eS z = e b b ~ a a with eigenvalues e(n& — n a )/2 and eigenstates 1^) = \ni, = v + 2k, n a = 
v + 2(M — k)). Here, the seniorities are the number of unpaired bosons v = 0,1, M = (j — v) is the total 




(32) 



ofEq. (14). 
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Figure 7: Real (a) and Imaginary (b) parts of the 15 pairons of the (k = 10) state of figure[6]as function of g(2j — 1). Vertical 
dashed lines indicate the g values where pairon collapses occur according to Eq.( |31| . The pairons close to t a = — t in the weak 
coupling limit remain trapped in the real interval [— t, t] , whereas the non trapped pairons expand in the complex plane after 
collapsing in the position Pp = 1 [horizontal gray line in panel (a)] of the effective charge Qb- The number of non trapped 
pairons for this 10-th excited state is 10, therefore collapses of 1 to 10 pairons occur at g values given, respectively, by g = 9jy 
with Np = 1, 10. The results corresponding to g^ _ 10 are indicated in the upper scale of the panels. For g < gjv D =io no 
more collapses occur, and the pairon set consist of 5 real pairons in interval [— t, t] and 5 complex conjugated pairon pairs. 



number of boson pairs, and k = 1, M. The positive (negative) parity sector corresponds to n a and rib 
even (odd), or in terms of the seniorities it corresponds to v = [y = 1). Independently of the parity, the 
ground state has M boson pairs occupying the a level. The excited sates are obtained by promoting boson 
pairs from the a to the b level. In this way the fc-th excited state for a given parity has M — k boson pairs 



occupying the a level, and k in the b level. From the wave function (19 1, it can be seen that in the fc-th 
excited state M — k pairons converge to t a = —t and k pairons to tb = +t. For finite but small g, it was 
shown in ref.|15j that the pairons e a can be approximated by 

e a « U ~ 5(1 + st 2 )r h (s = -l) (33) 

where U = t a = —t or ti — tb — +t, and ri are the positive roots of the Legendre polynomial V N , 1 ^ 2 (x), 
with N% (i = a, b) the number of e a pairons converging to ti for g = 0, i.e, N a = M — k and Nb = k, for 
the fc-th excited state. Hence, for small g, the entire set of states can be classified by the number of pairons 
distributed close to t a and close to tj,. For the fc-th excited state in the limit g negative and small, a group 
of M — k pairons is close and above t a = —t, implying they are in the interval [— t, t], which will turn be very 
relevant for their behavior in larger g. A second group of k pairons is close and above tb = +t, i.e. outside 
the interval [— t, t] . For small and positive g the situation is reversed. M — k pairons are outside the interval 
[— t, t] and close to t a = —t, whereas the remaining fc pairons are inside the interval and close to tb = +t. In 
Figure|6]we illustrate this behavior for the positive parity 10-th excited state of a system with M = j = 15 
pairons. In the figure, fc = 10 pairons are close to tb = +t and the remaining M — k — 5 sit close to t a = —t. 
Even if the perturbative result is not valid for large g, the behaviour of the pairons as a function of g is 
strongly dependent on their values at weak coupling. In particular on whether they are inside or outside 
the interval [— t, t]. The pairons inside the interval [— t, t] expand in the real axis but remain constraint to it 
for any value of g. On the contrary, the pairons outside the interval expand on the real axis moving away 
from the interval. For negative g this expansion takes place above the interval [— t, t] till the pairons start 



to collapse in the position Pjj = 1 of the effective charge Qo for values of g given by (31 ). Similar collapses 



occur for positive g for values given by (30), but here the collapses take place at the position Pc = — 1 of 
the effective charge Qc- 

For positive values of g, the fc-th excited state has M — fc pairon non trapped in the interval [— t,t]. 
Therefore, for the ground state (fc = 0) all the pairons will successively collapse at Pc — — 1 for increasing 
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g as it can be seen in figure [3] The first excited state has one pairon trapped in the interval [— t,t], while 
the other (M — 1) non trapped pairons successive collapse at Pc = — 1 for increasing values of g. The same 
reasoning extends for the rest of the excited states. For the most excited state, with all its pairons trapped 
in the interval [— t, t], no collapse occurs. 

The situation is somewhat reversed for negative values of g. The fc-th excited state has fc non trapped 
pairons. Therefore, the higher the excited state is, the more collapses will occur in the position of effective 



charge Qo at g values given by Eq. (31 1. For the grounds state (fc = 0), since all the pairons are trapped 
in the interval [— t, t] no collapse occurs, as it can be seen in figure [3] For the fc-th excited state successive 
collapses of 1 to fc pairons occur at g values given, respectively, by g% D= i, g% D=2 , ■■■,9N D =k- 

In Fig. [7] we illustrate this behavior for the same system and state of figure [6] (P — +, j — 15 and 10-th 
excited state), and for negative values of g. For small g, five pairons are located close and above t a = —t 
and fc = 10 pairons are close and above t\, — t. As g is increased the trapped pairons expand in the interval 
[— t, t], but remain constraint to it for any negative g. Whereas, the other ten pairons expand outside this 
interval, till they begin to collapse into the position (Po = 1) of the effective charge Qo at g values given 
by Eq.(31 1. These g values are indicated by vertical dashed lines in the panels of Fig. [7J Immediately after 



the collapse of an even number of pairons two complex conjugated pairons are created, and they expand in 
the complex plane until the next collapse. For the case illustrated in the figure, since the number of non 
trapped pairons is 10, the last collapse to Po — 1 takes place at g — g% D = w - From there on five complex 
conjugated pairon pairs expand in the complex plane for g < g^ =10 . 

A better insight on the pairon dynamics can be gained by plotting the pairon positions in the two 
dimensional complex plane. Fig. [8] shows the pairon positions in the complex plane of the complete set 
of positive parity states for a system of j = 30 and two different negative g. The fist one [panel (a)] is 
g = g% D= 22 where No — 22 pairons collapse in Po = 1- The second one [panel (b)] is an intermediate value 
of g between two collapses (5^0=23 < $ < 9n d =22)- The complex conjugated pairs of pairons of a given 
state are distributed in complex arcs around Po = 1. In panel (a), the radius of these arcs goes to zero for 
states with large enough number of non-trapped pairons, i.e. for those states with k = 22 to k = j = 30 non 
trapped pairons, corresponding to the 22-th to 30-th excited state. In both cases the outer arcs correspond 
to less excited states. 



Before closing this section, it is interesting to note that the condition (31) of No pairons converging to 



Pd = 1 in the hyperbolic region (s = —1) of phase diagram, defines hyperbolas in the 7^-7^ plane given by 

™* = G^+rLvJ ■ (34) 

For No = 1, the resulting hyperbola (7j,7:z = 1) is the same reported in [5] for the crossing of the ground 
states of positive and negative parities in the third quadrant. Moreover, in [31] it is argued that the rest 
of the hyperbolas (No — 2, M) define the points where there are crossings between excited states of the 
two parity sectors. For instance, for No — 2 the ground and first excited state of positive parity cross, 
respectively, those of the negative parity sector. In general, for arbitrary No the corresponding hyperbola 
defines the points where the first No states of positive parity cross, respectively the first No negative parity 
states. This result, already confirmed in [3TJI32] for the ground state, is numerically confirmed here for the 
ground and excited states in figure [9] The figure displays the absolute value difference {\Ep + — -Ep-|) 
between P+ and P— states for a system with j = 10 and t = 1/2 as a function of negative g (third quadrant 
in the phase diagram 7 X — 7^). The differences between the ground, first, second, third and fourth excited 
states of every parity sector are shown in logarithmic scale in order to make clear the crossings between 
positive and negative parity states. Every line is divided in continuous and dotted segments indicating if 
the difference Ep + — Ep- is negative or positive respectively. The points where this difference changes sign 
indicate a crossing bet wee n positive and negative parity states. The vertical dashed lines signal the points 



where the hyperbolas (34) are traversed, i.e. when g(2j — 1) = g c ND (2j — 1) = — 2j+i-2N D ' ^ s ^ can ^ e 
seen in the figure, for No = 1 [where g(2j — 1) = —1] the ground states of every parity sector cross, whereas 
for the leftmost vertical lines (A^ = 2,3,4,5), in addition to the ground states, more and more excited 
states cross. It is worth mentioning that the phase transition in the thermodynamic limit is expected at 
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Figure 8: Pairons in the complex plane for the complete set of positive parity states for a j = 30 system. Panel (a) corresponds 
to g = <?22 wnere collapses of 22 pairons are expected according to Eq. 1311. Panel (b) corresponds to a value of g where 
collapses are not expected. Horizontal dotted lines indicate the values t = —1/2, t = 1/2 and Pjj = 1. The complex pairons 
of each state accommodate in arcs around Pp = 1. Outer arcs correspond to lower energy states. In panel (a) the arcs of the 
k = 22 to k = j = 30 excited state collapse to the position = 1, with the outer arcs corresponding to lower energy excited 
states. 

g(2j - 1) = -t = -1/2 (see table |2j. We can appreciate in the figure a dramatic change in the difference 
between the energies of the positive and negative parity ground states around this value. However, the first 
crossing occurs at g(2j — 1) = —1. 

A preliminary view to the relation between collapses and crossings of different parity states indicates 
that the states participating in a given crossing do not have their pairons collapsing in Pp = +1. Contrarily, 
the pairon collapses of a given state prevent it from having a crossing. As result, the ground state having 
no pairon collapses, has crossings for all the values g^ . By contrast, the most excited state with pairon 
collapses in every value g c f fD , does not have crossings. For intermediate states, they present crossings in 
g = g^ D only if for this particular value, their pairons do not present collapses. This condition occurs if Njj 
is greater than the number of their non trapped pairons (k for the fc-th excited state). For instance, the 10-th 
P = + excited state of Figure [7] cross the 10-th excited P = — state only for 9pf D =u, 9%- D= x2> ■■■,9N r >=j- 
Further research in this relation is desirable to establish a deeper connection between both phenomena, 
crossings and pairing collapses. Likewise, the hyperbolic quadrants of the phase diagram have a region 
where avoiding crossings between states of the same parity take place. This region already identified in 
|33) by studying the density of states in the thermodynamic limit, appears in the hyperbolic case (s = — 1) 
when \g(2j — 1)| > 1/t. The relation between avoiding crossings and the pairon behavior in the SU(1, 1) 
RG models is out of scope of this contribution, and deserves more work for a complete clarification. 

6. Summary 

The Schwinger boson representation of the SU(2) algebra allows to connect the LMG model with the 
two-level SU(1,1) RG pairing models. We have exploited this relation to classify the entire parameter space 
of the LMG model in terms of the three RG families, the rational, the trigonometric and the hyperbolic. This 
classification sheds new light into the LMG phase diagram and its quantum phase transitions. Moreover, the 
electrostatic mapping of the trigonometric and hyperbolic models provides new insights into the structure 
of the different phases. We explored the LMG model from the perspective of the RG models, where the 
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Figure 9: Absolute value of the difference between the ground to the fourth excited state energies of positive and negative 
parity states for the hyperbolic LMG model with j = 10 and t = 1/2. Continuous segments indicate that the P = + state is 
lower in energy respect to the corresponding P = — state, whereas the dashed segments correspond to the opposite situation 
Ep— < Ep+. The vertical dashed lines indicate the values 9jv d (with Np> = 1, 5 signaled in the upper scale). 

eigenstates are completely determined by the spectral parameters (pairons) of a particular solution of the 
non-linear set of Richardson equations. We proposed a numerically robust method to solve the Richardson 
equations which generalize that of reference |28j . The method was proven to be suitable for obtaining the 
pairons of the complete set of eigenstates for moderate systems sizes. 

Using boson coherent states, we have re-derived the phase diagram of the LMG model and the character- 
istics of its different phase transitions. The second order phase transitions were interpreted in terms of the 
RG solution as a localization-delocalization of the ground state pairons, which takes place when the pairons 
concentrated around the value t a = —t at weak coupling, expand in the whole interval [t a ,tb\ — \—t,t\. On 
the other hand, the first order phase transition was related to the discontinuity of the pairons when the 
transition line is traversed. This discontinuity was related, in turn, with the crossings between states of the 
same parity. We have confirmed that the dynamics of pairons in the rational and trigonometric RG models 
take place entirely in the real axis. However, it was unexpected to find complex pairon solutions in the 
hyperbolic regions of the phase diagram. For the ground sate, complex values of the pairons are obtained 
after the collapses of an even number of pairons into the position Pc = — 1 of the effective charge Qc- It 
was numerically verified diagonalizing very large systems that no phase transition is associated with this 
singular pairon behavior, even for the particular case in which all the pairons collapse into Pc = — 1- For 
the latter case the ground state wave function has a particularly simple form which is equivalent to the 
Moore-Read state of the p x + ip y model [TH1 HH [20] ■ 

A complete classification of the excited states for the hyperbolic regions was given in terms of their 
pairon positions. For negative couplings it was found that the fc-th excited state is characterized by a set of 
(M — k) pairons trapped in the real interval [— t, t] while the other k pairons lie outside this interval. The 
dynamics of these non trapped pairons show collapses of Np pairons in the position Prj = 1 of the effective 
charge Qd, at g values given by g c N = —l/(2j + 1 — 2Ne>). This singular behavior of the excited state 
pairons was found to be connected to the crossings of Njj lowest positive parity states with Njj negative 
parity states, at exactly the same values cfc . As discussed in Ref . |22j , the pairon dynamics can help to 
identify significant physical phenomena. The relation between collapses and crossings is an example of this 
connection that deserves a deeper study. The collapses of pairons in the positions of the effective charges 
obtained in the hyperbolic boson RG models discussed here, is a feature also found in the p x + ip y fcrmion 
pairing realization of the hyperbolic SU(2) RG model (THJ, QjU [50] ■ In this latter model the collapses were 
related with another singular phenomenon: a third order phase transition. The insight gained in the study 
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of the LMG model, where the set of pairons of every state in the spectrum is easily accessible, can help to 
elucidate more intricate mechanisms in other integrable models where the numerical access to the pairon 
sets is more demanding. 
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Appendix A. Lame differential equation from the RG equations 

Here we show that the polynomial P(x) = Yl^ ( x ~ e ")' e a being the roots of the Richardson 

equations, satisfies the Lame differential equation. The demonstration is based on the simple identity 

1 < 1 1 ' (A.1) 



[x - F)(x -G) F-G\x-F x-G 

From the definition of P(x) it is straightforward to show that its derivative is P'(x) = X)^Li T\p^ a ( x ~ e a)- 
Therefore, we have the identity A(x) = P'(x)/P(x) = J2 a x -e ■ ^ derivative of A(x) yields 

1 ' P(x) y ' ^ {x-e a y 
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From the previous result and the identity ( |A.1[ ), we obtain 

P"{x) l ^l^l 



E x - e D E 



e/3 



Now, if e a are the roots of the Richardson equations, the second sum in the right hand can be substituted 
according using Eq. ( 22 1 



P"{x) 
~P{x) 



E^E 



1 



1 



where we have, again, used the identity ( A.l ). From the definition of A(x) the previous equation is equivalent 
to 



P"{x) 
~P(x) 



E 



2p k \ P'(x) 
x-rjk P(x) 



E 



2p k A(i] k ) 



= 0. 



Finally, multiplying the previous expression by P(x) and A(x) = Yik( x ~ we obtain the Lame differential 
equation (24), 



A(x)P"(x) + 



k=l 



2p k 

x - r) k 



P , (x)-[^PkHVk)'[[^ 



l^k 



P(x) = 0. 




Appendix B. Condition to have Nc pairons converging in the position Pc — —1 of the 
effective charge Qc in the hyperbolic LMG model 

Let us assume that N c pairons converge to Pc = — 1- Therefore we expand them as e a = — 1 + Sz a , with 
a = 1,...,Nq, and 5 an infinitesimal parameter. In this limit the Richardson equations corresponding to 
a = 1, ...,Nc separate in a term proportional to 1/5 and terms of order O(S ). Both terms have to cancel 
independently. The term proportional to 1/8 reads 



with a = 1, ...,N C 



Therefore, in order to have Nc pairons converging to Pc — — 1 the previous set of equations has to have a 
solution. Following the lines of | Appendix A[ let us suppose a polynomial whose roots give the solution of 
the previous set of Nc equations, Pn c ( z ) — Yia=i( z ~~ 7 -a)- If the set z a solves the system of equations, it 
can be shown that the previous polynomial is a solution of the differential equation 

zP^ c (z) + 2P' Nc (z)Q c = 0. 

The general solution is Pn c ( z ) — D \-2Qc — ^ w ^ n ^ arL0 - ^ integration constants. By comparing this 
solution with the initial assumption [Pn c = T\H=i( z ~ z a)], we obtain the conditions D = (1 — 2Qc) and 
— 2Qc + 1 — Nc- From latter condition and the definition (21 1 of Qc we finally obtain that the necessary 
condition to have Nc pairons converging into Pc = — 1 is 



9 = 9n c 



2j + 1 - 2N C 

Note that this condition applies to any eigenvector, ground or excited state. Following a similar reasoning it 
can be shown that the condition to have Nr> pairons converging to Po = 1 is —2Qn + 1 = No, from where 
the condition 



9 = 9 N D 



1 



2j + l- 2N D 

follows. It is important to note that the latter condition is satisfied for negative g, whereas the former is for 
positive g. 
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